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Dominant QCD Backgrounds in Higgs Boson Analyses at the LHC: 
A Study of pp — > ti + 2 jets at Next-To-Leading Order 
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We report the results of a next-to-leading order simulation of top quark pair production in asso- 
ciation with two jets. With our inclusive cuts, we show that the corrections with respect to leading 
order are negative and small, reaching 11%. The error obtained by scale variation is of the same 
order. Additionally, we reproduce the result of a previous study of top quark pair production in 
association with a single jet. 
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Exploiting the relatively clean experimental signal and 
excellent theoretical understanding of the W + W~ Higgs 
boson decay channel, the Tevatron collaborations con- 
tinue excluding Higgs boson masses, Mh, in the vicinity 
of twice the W boson mass with ever growing confidence 
levels [lj. As a consequence, studies assuming a lighter 
Higgs boson are more timely than ever. Clearly, it is nec- 
essary to devise suitable scenarios that would allow for 
discovery and measurement of the basic properties of the 
scalar at the Large Hadron Collider (LHC). Among the 
many available scenarios, one is especially disputed. It 
has been argued that if My < 135 GeV, the production 
and decay chain pp — > tiH* — > ttbb should provide an 
excellent opportunity owing to the large top and bottom 
quark Yukawa couplings. The question of discrimination 
between the large backgrounds and signal was expected 
to be solved by the sharp Breit-Wigner peak of the Higgs 
boson in the invariant mass distribution of the b quark 
jets. Unfortunately, application of realistic selection cuts 
and acceptances demonstrated a smearing of the reso- 
nance far beyond what would be expected from initial 
state radiation @, |||. At this point it seems, therefore, 
that a very precise knowledge of the backgrounds is nec- 
essary, if the channel is to be of any usefulness 23 . 

A close scrutiny of the backgrounds shows [| , Q that 
the most relevant are the direct production of the final 
state ttbb (irreducible background) and the production 
of a top quark pair in association with two jets, ttjj 
(reducible background). The latter needs to be taken 
into account due to the finite efficiency in identifying fa- 
quarks in jets (b-tagging) . We depict example diagrams 
contributing to the signal and the two backgrounds in 
Fig.[TJ It is crucial that although ttjj has a cross section, 
which is larger than that of ttbb byabout two orders of 
magnitude, in the actual setup of [3j, both backgrounds 
turn out to give very similar contributions. Recent stud- 
ies 0-0] have demonstrated that the next-to-leading or- 
der QCD corrections to ttbb are very large, undermining 
even further the feasibility of actual analyses in the chan- 
nel at hand. A proposed solution Q involves the use of a 



FIG. 1: Example diagrams contributing to the signal tiH* — > 
ttbb, and the irreducible background with the same final state, 
as well as the reducible background with two jets. Thick blue 
lines correspond to top quarks, the wiggly to gluons as usual. 



dynamical renormalization/factorization scale in simula- 
tions and, more importantly, the imposition of a jet veto 
on the third jet (upper bound on the allowed transverse 
momentum, pr). The purpose of this letter is to deter- 
mine, whether large corrections affect the ttjj channel. 

In principle, the concept of a cross section of a pro- 
cess involving massless partons at leading order (LO), 
as is the case of ttjj, is undefined due to potential soft 
and collinear divergences and requires the specification 
of separation cuts. This complicates somewhat our prob- 
lem, since any conclusions we will draw will pertain to the 
particular setup we will have used. Therefore, in order 
to retain some level of generality, we will not attempt to 
reproduce the exact conditions of the Higgs boson anal- 
ysis of Q. On the contrary, we will mimic some of the 
general assumptions, but impose more inclusive cuts. 

To be specific we consider proton-proton collisions at 
the LHC with a center of mass energy of y/s = 14 TeV. 
We set the mass of the top quark to be m t = 172.6 
GeV and leave it on-shell with unrestricted kinemat- 
ics. The jets are defined by at most two partons using 
the kx algorithm of d, 0] with a separation AR = 0.8, 
where AR = J(m - y 2 ) 2 + (fa - c/> 2 ) 2 , V < = 1/2 HE, - 

the azimuthal 



v/(2/i-y 2 ) 2 + (0i-02) 2 , 

Pi, z )/(Ei +Pi,z) being the rapidity and 
angle of parton i. Moreover, the recombination is only 
performed if both partons satisfy \yi\ < 5 (approximate 
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detector bounds). We note that the kr algorithm spec- 
ifies not only which partons are combined into jets, but 
also the momentum of the resulting jets. In our simula- 
tion, we assume that the four- momenta of the partons are 
added. We will comment on the importance of this point, 
when we compare with the tij calculation of [13, LUJ. We 
further assume that the jets are separated by AR = 1 
and have |j/j t| < 4.5. Their transverse momentum is re- 
quired to be larger than 50 GeV. It is mostly in AR and 
PT.min that the present setup is different from that used 
in the case of tibb [B|, Q . Here, we use a higher px cut 



similarly to the ttj case of The additional separa- 



tion in AR is only used to demonstrate the flexibility of 
our tools, but is believed to bear no impact on the final 
conclusions. Notice, that our jets are allowed to contain 
b quarks, and therefore also the final state tibb. This is 
irrelevant, since we know that this contribution is tiny 
and can be simply subtracted from our results (the sum- 
mation is incoherent). Finally, we note that the third jet, 
which stems from real radiation, is not restricted. Nev- 
ertheless, we will also study the impact of a jet veto with 
p T = 50 GeV. 

The production of jets in hadronic collisions can be 
decomposed into many processes at the parton level. In 
Tab. [H we summarize the LO contributions of different 
subclasses. These were obtained with the kinematics 
specified above using the CTEQ6L1 LO parton distribu- 
tion functions 12|, |l3[ and the LO running of the strong 
coupling constant up to the scale (common for renor- 
malization and factorization) fj,Q — m t- There are two 
interesting conclusions to be drawn here. First, it comes 
as a surprise that the mixed channel qg is more impor- 
tant than gg. The proportions between the two depend, 
however, on the px cut. If pT.mm = 20 GeV, the situ- 
ation is reversed. The reason is that the lower the cut, 
the larger the soft gluon enhancement. If there are more 
gluons in the final state the enhancement will be higher, 
and therefore at some point gg — > tigg must be larger 
than qg — » ttqg. The second, more important, point is 
that the channels related to the tibb final states, in par- 
ticular gg — > tiqq arc almost negligible compared to the 
two dominant (we referred to this before). This implies 
that whatever the result on the size of the corrections in 
the tibb case, tijj requires a separate study. 

Before we give our results for the next-to-leading order 
(NLO) corrections, we are compelled to present the com- 
putational framework used for the simulations. Similarly 
to our previous publication Q , we have used the Helac- 
Phegas framework lil-tlrj] and in particular Helac-IL 
17| for the evaluation of the virtual corrections. This tool 
uses the OPP method [H and CutTools [HHUI for the 
reduction of tensor integrals, as well as OneLOop for 
numerical values of scalar integrals. The practical tech- 
niques involve re-weighting of events and sampling over 
polarization and color, for details see [6]. The real ra- 
diation corrections were evaluated with Helac-Dipoles 



Process 


„LO LL] 


Contribution 


PP -> tijj 


120.17(8) 


100 % 


qg ->■ tiqg 


56.59(5) 


47.1 % 


gg -> tigg 


52.70(6) 


43.8 % 


qq' — !> tiqq , qq — ¥ tiq'q' 


7.475(8) 


6.2 % 


gg ->■ tiqq 


1.981(3) 


1.6 % 


qq ->■ tigg 


1.429(1) 


1.2 % 



TABLE I: The LO cross section for pp — > tijj production at 
the LHC. The individual contributions of the various partonic 
channels are also presented separately. Both q and q' span all 
quarks and anti-quarks. 



^22j . which is an implementation of the Catani-Seymour 
subtraction formalism [3, [24| . In order to check our re- 
sults, we have explored the independence of the results 
on the unphysical cutoff in the dipole subtraction phase 



space, see 22] and references therein for details. We have 
also verified the cancellation of divergences between the 
real and virtual corrections. Finally, the numerical pre- 
cision of the latter was assured by using gauge invariance 
tests and use of quadruple precision. We note that the 
only new virtual amplitudes are those involving a top 
quark pair and four gluons. These were presented for the 
first time in 17[ and, due to their notorious complexity, 
still await an independent check by other groups. 

For the evaluation of the NLO corrections, we have 
used the CTEQ6M parton distribution functions with 
NLO running of the strong coupling constant. At the 
central scale /iq — m t , we obtain 



(7 



NLO 

pp^ttjj+X 



(106.94 ±0.17) pb , 



where the error comes from Monte Carlo integration. 
Compared to the LO result from Tab. U this represents 
a negative shift of 11%, and allows us to conclude that 
the corrections to this process are small. 

The scale dependence of the corrections is illustrated 
in Fig. [2] At first, we observe a dramatic reduction of the 
scale uncertainty while going from LO to NLO. Varying 
the scale up and down by a factor 2 changes the cross 
section by +72% and -39% in the LO case, while in the 
NLO case we have obtained a variation of -13% and - 
12%. Second, the central scale that we have chosen is 
very close to the point of minimal corrections and slightly 
above the point of maximum of the NLO cross section. 
Indeed, both /i = 1/2^,q and /i = 2^o give smaller values. 

Taking into account the above dependence on the scale 
choice, it is to be expected that adding a jet veto will only 
worsen the result. In order to make this more transpar- 
ent, it is best to consider the difference between the cross 
section without and with the jet veto. This difference is 
given by a LO calculation, since it requires the existence 
of three separated jets with a lower cut on their respec- 
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FIG. 2: Scale dependence of the total cross section for pp — > 
tijj + X at the LHC with /j,r = /if = £ ■ fJ,o where /io = TUt- 
The blue dotted curve corresponds to the LO, the red solid to 
the NLO result whereas the green dashed to the NLO result 
with a jet veto of 50 GeV. 



tive transverse momenta. Clearly, not only will it be neg- 
ative, but it will also have to grow more negative with 
diminishing scale (since this is almost entirely governed 
by the behavior of the strong coupling constant alone). 
At this point it is difficult to decide whether a jet veto 
will or will not be necessary for the complete Higgs bo- 
son analysis. What we can give is the lower bound on the 
corrections, assuming that a jet veto of less than px — 50 
GeV is likely to endanger the stability of the perturba- 
tive expansion. The total cross section with a jet veto of 
50 GeV is 



NLO / 



> T ,x < 50 GeV) = (76.58 ± 0.17) pb , 



with a scale variation of -54% and -0.3%, see Fig. [5] The 
plots show that choosing a higher scale in the case of 
the jet veto, would lead to a result with virtually no 
scale dependence. This should be considered as severely 
underestimating the error. 

While the size of the corrections to the total cross sec- 
tion is certainly interesting, it is crucial to study the cor- 
rections to the distributions. The most important for us 
is the invariant mass of the two tagging (highest pr) jets, 
since this is the observable entering Higgs boson studies. 
We plot the LO and NLO results in Fig.H While we no- 
tice a long tail, we keep the dependence only in a modest 
range up to 400 GeV due to our phenomenological moti- 
vation. The distribution starts above about 45 GeV due 
to the AR and pr cuts, and shows tiny corrections up to 
at least 200 GeV, which means that the size of the cor- 
rections to the cross section is transmitted to the most 
relevant distribution. 

Of course, there are observables showing much larger 
effects. The classic example is the transverse jet mo- 
mentum distribution at high pr- We illustrate the 
phenomenon in Figs. U and [5J which demonstrate the 
strongly altered shapes in the cases of the hardest and 



> 

0) 

O 

\ 

ft 



lO- 



T3 



10 



-2 




100 



200 



300 



400 



mjj [GeV] 



FIG. 

and the second hardest jet for pp —¥ tijj + X at the LHC. The 
red solid line refers to the NLO result while the blue dotted 
line to the LO one. 
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FIG. 4: Distribution in the transverse momentum pTj of the 
1st hardest jet for pp — >• tijj + X at the LHC. The red solid 
line refers to the NLO result while the blue dotted line to the 
LO one. 



second hardest jets. It is well known that this kind of 
corrections can only be correctly described by higher or- 
der calculations. On the other hand, the behavior at low 
Pt is certainly further altered by soft-collinear emissions, 
which are best simulated by parton showers. With our 
lower cut of PT,min — 50 GeV, we expect to be mostly 
in the safe range, where fixed order perturbation theory 
does not break down. 

While the above comments conclude our analysis of 
ttjji we would like to make a few statements about the 
process with only one jet, ttj. While preparing our cal- 
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FIG. 5: Distribution in the transverse momentum pr } of the 
2nd hardest jet for pp — > tijj + X at the LHC. The red solid 
line refers to the NLO result while the blue dotted line to the 
LO one. 



culation, we have made a comparison with the results of 
11]. Using exactly the same setup as that work, we were 
able to obtain the following value for the cross section 



a 



NLO 

pp—yttj+X 



(376.6 ±0.6) pb 



demonstrating perfect agreement ([ll( quotes (376.2 ± 
0.6) pb). We note that besides a different value of the 
top quark mass (m t = 174 GeV) and Ai? = 1 in the jet 
algorithm, the authors used a kx algorithm, where the 
momenta of the partons being combined into a jet are 
not simply added. Instead, massless jets are constructed 
using transverse momenta (which are trivially added) and 
rapidities. The version that we use would lead to a lower 
value of the cross section, since the transverse momentum 
of the sum of momenta is never larger than the sum of the 
transverse momenta, and the only relevant cut is applied 
to the jet pr- Our value would instead be (372.2±0.6) pb. 
On the other hand, with all parameters and cuts as given 
in the introduction the final result is (376.1 ± 0.7) pb. 

Let us conclude by pointing out that while we demon- 
strated that the corrections in tijj are small for one 
setup, comparison to the behavior of the corrections for 
different setups in the ttj case 



11 



provides a viable argu- 
ment that our conclusions will remain true for other input 
parameters. Nevertheless, we are planning to present a 
much wider study involving in particular a variation of 
the center of mass energy, cone size in the jet algorithm, 
transverse momentum cuts and jet vetoes. 

The simulations presented in this work have been per- 
formed at the DESY Zeuthen Grid Engine computer clus- 
ter. The work of M.C. was supported by the Heisenberg 
Programme of the Deutsche Forschungsgemeinschaft. 
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